Conical intersections in an ultracold gas 
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We find that energy surfaces of more than two atoms or molecules interacting via dipole-dipole po- 
tentials generically possess conical intersections (CIs). Typically only few atoms participate strongly 
in such an intersection. For the fundamental case, a circular trimer, we show how the CI affects 
adiabatic excitation transport via electronic decoherence or geometric phase interference. These phe- 
nomena may be experimentally accessible if the trimer is realized by light alkali atoms in a ring trap, 
whose dipole-dipole interactions are induced by off-resonant dressing with Rydberg states. Such a 
setup promises a direct probe of the full many-body density dynamics near a conical intersection. 
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Introduction: Conical intersections (CIs) of electronic en- 
ergy surfaces are a generic feature of large molecules 
[H, The intersections provide fast intra- molecular 
transition channels between electronic states and can 
thus affect the outcome of (photo-) chemical processes. 
Radiation-less de-excitation of large bio-molecules pro- 
ceeds through these channels and yields enhanced pho- 
tostability that might have been crucial for the devel- 
opment of life on earth Modern techniques allow 
the theoretical investigation of quantum wave-packet dy- 
namics at conical intersections of large molecules in quite 
some detail 0, 0] • However experiments usually monitor 
such dynamics indirectly, e.g., through reactant fractions 
or fluorescence spectra. 

The direct observation of many-body densities could 
be realized near conical intersections in ultracold gases. 
As we demonstrate, these are ubiquitous between energy 
surfaces of inter- atomic (or molecular) dipole-dipole in- 
teractions [8]. Such interactions are responsible for exci- 
tation transport in photosynthetic light harvesting units 
molecular aggregates [loj or cold Rydberg gases 11 1. 
We find that typical ensembles have sub-units of three 
to seven particles that are mainly responsible for the CI. 
We investigate the fundamental unit, a ring-trimer, in 
more detail. In the ring geometry, two adiabatic energy 
surfaces cross at the equilateral triangle configuration, 
forming a conical intersection. While a direct crossing of 
the CI splits wave packets among two electronic surfaces 
and thus de-coheres the reduced electronic state, circum- 
navigation of the CI results in interference patterns with 
clear signatures of the geometric phase. 

We suggest a realization of the ring-trimer using 
trapped atoms that obtain dipole-dipole interactions 
through off-resonant Rydberg dressing [12|. Through 



spatially resolved detection of single atoms one can re- 
veal the full many-body dynamics around the CI. 
Conical intersections in dipole-dipole gases: Consider a 
sample of N atoms or molecules of mass M whose po- 
sitions r n are grouped into the 3N component vector 
R = {r n }. We assume that each particle can be ei- 
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FIG. 1: (color online) Conical intersections in ensembles of 
dipole-dipole interacting atoms or molecules, (a) Randomly 
placed particles (*) and location of a selected conical inter- 
section (red •) in their vicinity (grey shade, R ma x — 0.2d m ax, 
see text), (b) Potential energy surfaces Uk(Av\ r + AiuW) 
near the CI in (a) using a = 1. The modes V (W) shown 
as red (black) arrows in (a) span the CPs branching plane, 
see Eq. (|3]). (c) Mean connectivity iVconn (solid line), as ex- 
plained in the text, and its standard deviation (dotted lines), 
(d) Histogram of participation numbers for Rmax = 0.35d m ax. 



ther in some electronic ground state \g) or some ex- 
cited state \h). The single excitation Hilbert space 
for the electronic degree of freedom is then spanned by 
| 7r n ) = \g • • • h - ■ ■ g), with the n'th particle excited. 
We incorporate particle motion and electronic excitation 
transfer by dipole-dipole interactions [13| into the Hamil- 
tonian 
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where F nm is the dipole-dipole interaction between the 
particles. Besides Rydberg atoms 14j, LL5j, this model 
can also describe molecular aggregates [10j or Rydberg- 
dressed ground state atoms |l2j. We ignore the angu- 
lar dependence of the interaction and specify V nm — 
—fi 2 /\R nrn \ 3 with Rnrn = |r n — r m |. Atomic units are 
implied unless otherwise indicated. 

To analyze the Hamiltonian (pQ), consider a Born- 
Oppenheimer separation of dynamics [l6|. The 
eigenstates of the electronic Hamiltonian fulfill 
# el (R)|(^-(R)) = ^-(R)I^-(R)). The energies 
Uj (R) then define TV adiabatic potential energy surfaces 
for particle motion. A particle stays on a definite 
surface, if non- adiabatic couplings between the surfaces 
can be ignored. These non-adiabatic couplings become 
large at a conical intersection. 

To assess the relevance of conical intersections in an 
ultracold gas, we consider an ensemble of N atoms dis- 
tributed randomly in a cube of side length L, with coor- 
dinates R* and an isotropic quantum or thermal position 
uncertainty of about or* . After determination of R* , we 
seek Cl-locations Hci within a distance R ma x = &r* [I3 
by numerical minimization of the smallest separation be- 
tween any two surfaces U\ and Uj. We define N conn as 
the number of adjacent energy surfaces between which 
we find a CI under these constraints. Its average is 
shown in Fig. [T] (c) as a function of Rmax/dmax, where 
d m ax = L/(2ttN) 1 ^ 3 is the peak of the nearest neigh- 
bor distance distribution. We averaged over 50 random 
atomic ensembles and 200 different minimization seeds. 
Already for Rmax ~ 0.2dm ax ~ 0.06L at least 3 CIs are 
accessible on average. An ensemble with spread of posi- 
tion or* ~ 0.2d max , would necessarily access these CIs. 

To characterize CIs between surfaces i and j in more 
detail, consider their branching plane spanned by [l| 

2vj? = (tpi \V Tn Hel\<Pi > + <¥>j \V rn Hel\<Pj ), 
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Let = {vjf }, W y ' = {w^} be 3N component vec- 
tors and V^' = V^'/IV^'I, W^' = /\W^\. The 
motion of atom m has a large effect on the energy gap 
between surfaces i and j, if the |v m | 2 or |w m | 2 compo- 
nents of the branching vectors are large. Motion in the 
seam-space orthogonal to V ZJ , W ZJ has almost no effect, 
thus small |v m | 2 or |w m | 2 indicate that the atom does 
not participate strongly in the intersection. This moti- 
vates the use of a participation number (see e.g. [Uj) 
N'J = (VE„|v^| 4 + l/EJw^'| 4 )/2, which roughly 
gives the number of atoms involved in CI dynamics. The 



three to seven. Changes in density p = N/L 3 only affect 
the overall energy scale of Fig. [T] (b) but not the results 
in Fig.[U(c-d) if distances are expressed in units of <i m ax- 

We conclude that conical intersections are generic in 
dipole-dipole interacting aggregates with moving con- 
stituents. This holds also for tightly bounded motion, 
such as small vibrations around fixed locations. Since a 
small number of atoms control the CIs we will study the 
fundamental case N = 3 in the following. 
Ring trimer: We consider three particles confined one- 
dimensionally on a ring, with coordinates as in Fig. [2] On 
the ring the position of atom n is specified by its 2D polar 
angle 6 n . Then the dipole-dipole interaction is V nm = 
— ji 2 / (R^/2(l — cos[# n — m })) 3 . The centre-of-mass type 
angle 0cm = ^Z n ^n/3 decouples, and the dynamics of 
interest is fully described by the two relative angles 9 12 = 
6> 2 -6>i and 6> 23 = 3 - 2 . 

In Fig. [2] we show the three potential energy surfaces 
obtained by diagonalizing H G \. The top (red) and bot- 




FIG. 2: (color online) Ring trimer. (a) Three atoms con- 
fined to a ring with radius R in the x-y plane, (b) Born- 
Oppenheimer surfaces for the electronic Hamiltonian (pQ). 
From top to bottom Urep, ^mid U a tt> For better visibility 
we shifted the repulsive (attractive) surface up (down) by 
AU = 0.37 MHz, and magnified the middle surface by a fac- 
tor fm = 5. At the marked CI location, Urep and U m [^ would 
touch without the shifts. 

torn (blue) surfaces are globally repulsive or attractive, 
respectively. Each of them is associated with a different 
delocalized excitation [14[. Energetically between these, 
we find a middle surface that touches the repulsive one 
in a conical intersection at #12 = #23 = 27r/3. In the 
following we show how this intersection affects adiabatic 
excitation transport. 

The total quantum state is | ^(R) ) = 
Z)n=i 0n(R)| TTn ), where n (R) is the wave func- 
tion of particle motion in the electronic state 1 7r n ) . 
From the Hamiltonian (pQ) one obtains the equations 
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A/p-histogram in Fig. [T] (d) shows that these are usually For N = 3 atoms we use #12 and #23 and write | ^ ) = 
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J2n=i <f>n(0i2, ^23)| ), which reduced Eq. (gj) to 
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We also introduce the adiabatic form of the wavefunc- 
tion I*) = E?=i&(*i2, 023)1 ip j (6 12 ,6 23 )). Finally, let 
us define the population of the n'th diabatic (adiabatic) 
state as p n = J dH\(j) n \ 2 (pj = j dR|^| 2 ). For adia- 
batic states and surfaces, we will use subscripts rep = 1, 
mid = 2, att = 3 as in Fig. [2j The setup described 
above allows one to observe how CIs affect many-body 
dynamics. We will discuss two paradigmatic situations. 
Electronic de- coherence: A robust consequence of the CI 
is to decohere reduced electronic state of the the trimer. 
Consider the case where initially the system is on the re- 
pulsive adiabatic surface and the particles form an isosce- 
les triangle as skeched in Fig. [2] (a). The initial dis- 
tance between particles 1 and 2 is R12 = do < <i*, where 
d* = V3R is the separation at the conical intersection. 
Each atom has initially a Gaussian angular distribution 
with width <jo <C 2tt. The detailed construction of the 
initial state is described in [15j. From this initial state, 



the repulsion drives the configuration over the CI. We 
use R 9.8/im, /i 180 a.u. and M 11000. These 
parameters correspond to laser dressed Lithium, as will 
be explained later. 

The close proximity pair of particles initially accel- 
erates adiabatically, as evident from the constant sur- 
face populations in Fig. [3] (b). However, as the many- 
body wave function passes the conical intersection, more 
than 60% of the total population is transferred from 
the repulsive to the middle surface. Fig. [3] (c) shows 
an accompanying drop in purity of the reduced elec- 
tronic density matrix a = Z) n ,m a nm| 7r n )( 7r m |, with 
o~nm — J d N H 0*(R)0 m (R). The operator a describes 
the electronic state, disregarding knowledge of particle 
positions. The observed decoherence of the electronic 
state is a consequence of the spatially disjunct splitting 
of the wave-packet on two surfaces and ensuing entangle- 
ment between the particle position and aggregate elec- 
tronic state. It will occur generically for all dynamics 
close to the CI. 

Geometrical phase: A more subtle but even more funda- 
mental consequence of conical intersections is the geomet- 
ric (Longuet-Higgins-Berry) phase picked up by wave- 
packets that are adiabatically transported in a closed- 
loop around them [l9|, [lo|. In our case, wave-packets 
moving half a circle around the CI on the repulsive sur- 
face with opposite sense of rotation meet with a relative 
phase of tt. We initially superimpose two different mo- 
mentum components, which yield the classical trajecto- 
ries shown in Fig. 3] (a) as solid white lines. The resulting 




FIG. 3: (color online) (a) Total atomic density, given 
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E n =J d N ~ 1 6 {m} \(f) n (6)\ 2 , where / d N ~ 1 6 {m} denotes inte- 
gration over all but the ra'th particle coordinate. Note that 
the #-axis is periodic, hence the highest and lowest parti- 
cle are near neighbors and initially repel. The configuration 
passes the CI location when all three inter-particle distances 
are equal, around t — 0.17 ms. (b) Adiabatic populations 
Prep (solid, black), p m id (solid, red), electronic populations 
Pi — P3 (dashed, blue), P2 (dashed, red) and total popula- 
tion (dotted, black), (c) Purity P = Tr[a 2 ] of the reduced 
electronic density matrix, described in the text. 



quantum interference pattern after encircling the CI is 
shown in Fig. |4] (b) for R = 13/im. The wave-packet is 




FIG. 4: (color online) Intersection circumnavigation from ini- 
tial momentum superposition. The resulting interference pat- 
tern manifests the geometrical phase, (a) Shown is the den- 
sity |0re P | 2 at times ti = 0, t 2 = 0.24ms, £3 = 0.48ms. The 
white lines are classical trajectories on the repulsive surface. 
Segregation of different time samples is provided by the white 
dashed ellipses, (b) Magnification of the area within the violet 
square in (a). A 7r-shift of the pattern due to the geometrical 
phase is visible between the part of the wave packet that has 
circumnavigated the CI (<g>) and the part that has not. Black 
lines in the background are iso-contours of U rep - 



chosen to contain components that will not actually have 
enclosed the CI upon interference. These form the lower 
left part of the pattern in Fig. 2] (b), with an antinode 
on the diagonal. In contrast, the top right of the pattern 
resulting from motion enclosing the CI has a node on the 
diagonal. This manifestation of Berry's phase shift can 
be experimentally accessed due to the control and detec- 
tion afforded by ultracold atoms as we will show in the 
next section. The interference pattern could then be used 
as a sensitive probe of the many-body energy landscape, 
as for example the presence of a fourth perturber atom 
would move the CI out of the #12, #23 plane and hence 



yield geometric phases differing from tt [l|. 
Laboratory realisation: We propose a realisation of the 
ring- aggregate model using trapped alkali atoms. To in- 
crease flexibility in the choice of parameters, we assume 
that long-range interactions are induced by dressing with 
Rydberg states [Hj], which allows the atoms to be pre- 
dominantly in their ground-state. For this we require two 
non- interacting ground (or meta-stable) states \g), \h), 
each of which is off-resonantly coupled to a Rydberg level 
\g) | s ), \h) \p) by two separate transitions, both 
with Rabi- frequency Q and detuning A. The detuning 
parameter a = Q/2A controls the population transferred 
to the Rydberg levels | s), \p), of order ~ a 2 , and the 
interaction-strength for the ground state atoms (~ a 4 ). 

Let the Rydberg levels \s), \p) have a principal quan- 
tum number v. We parameterize the dipole-dipole in- 
teraction strength \i and life-time r of the dressed Ryd- 
berg states as \i = \i^v 2 o? and r = tqv 3 / a 2 respectively, 
where fio and To can be found from a reference Rydberg 
level. We will explicitly consider 7 Li with M = 11000, 
Ho = 0.8, To = 3 x 10 T [21] (all in atomic units). 

The simulations presented above correspond to 7 Li 
atoms confined on a ring of radius R = 9.8/im, acquir- 
ing long-range interactions through laser dressing with 
a = 0.15 via the Rydberg level v = 100. The parameters 
have been chosen to meet the following conditions: 

1. The dipole-dipole interaction energy does not exceed 
realistic transverse oscillator spacings uj± of the ring, en- 
abling one- dimensional dynamics. The effects presented 
can occur close to the CI where the atoms are sepa- 
rated by do « d*, their separation at the CI. Their 
dressed interaction energy is then V(d*) = p 2 /dl = 

/(V3R) 3 = 40 kHz. This is within reach of very 
tight trapping potentials. 

2. The bare dipole-dipole shift without dressing is smaller 
than the dressing laser detuning, at the closest approach 
of the atoms. Otherwise resonances would yield undesir- 
ably large excited state populations. We have Vbare(^*) = 
V(d*)/a 4 = 80MHz_^ while a detuning of A ^ 150 is fea- 
sible for z/ = ioo da. 

3. The motion is adiabatic, except close to the conical 
intersection, as seen in Fig. [3] 

4- The life-time of the dressed three- atom system (r = 
rov 3 /a 2 /N = 10ms) exceeds the duration of the motion. 
A systematic search for regions in parameter space where 
realistic ring-trimers exist is non-trivial and will be re- 
ported elsewhere. It shows that light alkali atoms such 
as Li are favorable and that dressing is crucial. Without 
it, the above conditions can only be fulfilled for unrealis- 
tically large principal quantum numbers. 

In conclusion, we have shown that external conical in- 
tersections in dipole-dipole interacting gases of ultracold 
(Rydberg) atoms or molecules are common, as in the in- 
ternal energy landscape of large molecules. We have pre- 
sented a setup which holds promise for direct laboratory 
studies of the many-body wave function dynamics near 



or across a conical intersection. This can be achieved in 
a circular aggregate of Rydberg-dressed alkali atoms con- 
fined in a ID ring geometry, for which we have worked 
out paradigmatic consequences of the CI such as elec- 
tronic decoherence and nuclear interference due to the 
geometric phase. Conical intersections in laboratory co- 
ordinates can also be designed through external fields for 
23| , but lack by construction the 



ultra-cold molecules 22 



many-body interacting character of the CIs studied here. 

Similar effects as in our examples would occur also if 
the dipole-dipole interacting particles were not free on a 
ring, but confined in separate potential wells, as long as 
vibrations within each well allow many-body wave packet 
excursions over the CI. Hence, our results pertain in par- 
ticular to molecular aggregates such as light harvesting 
units. 
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